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Abstract 

The small- and large-sphere limits of the quasi-local energy recently proposed by Liu and Yau 
are carefully examined. It is shown that in the small-sphere limit, the non-vacuum limit of the 
Liu-Yau quasi-local energy approaches the expected value -^r 3 X(eo, eo). Here, % is the energy- 
stress tensor of matter, eo 6 T p M is unit time-like and future-directed at the point p located at 
the center of the small sphere of radius r in the limit r — > 0. In vacuum, however, the limiting 
value of the Liu-Yau quasi-local energy contains the desired limit |^Q5(eo, eo, eo, eo), where <B is 
the Bel-Robinson tensor, as well as an extra term. In the large-sphere limit at null infinity, for 
isolated gravitational sources, the Liu-Yau quasi-local energy is shown to recover the Bondi mass 
and Bondi news flux, in space-times that are asymptotically empty and flat at null infinity. The 
physical validity of the Liu-Yau model in view of these results is discussed. 
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I. INTRODUCTION 



Among a handful of unsettled puzzles at the foundation of Einstein's general theory of 
relativity, the very basic notion of energy- momentum seems to be of everlasting interest. 
Despite the triumph of the proof of the positivity of the total gravitational energy at both 
spatial and null infinity [31|, there is a lack of a well-defined notion of the local gravitational 
energy-momentum density. In fact, the equivalence principle, or, the existence of the normal 
co-ordinate system, prohibits any non-trivial point-wise localizable density. Consequently, 
only quasi-local quantities are meaningful. Although various proposals have been put 
forward (see, for example, (32] for a fairly complete and up-to-date review.), it should be 
noted that the study of quasi-local quantities is still rather premature in the sense that no 
truly axiomatic framework has been distilled from physics. A generally accepted strategy of 
studying quasi-local quantities is to devise such quantities and to check that they recover, 
in certain limiting situations, known properties. To be more specific, take for example the 
quasi- local energy-momentum. It has been proposed 8j, 32|] that it satisfy the following 
empirical criteria: 



(CI) Causality: 

Quasi-local energy-momentum is future-directed and non-space-like, provided that matter, 
if any, satisfies the dominant energy condition in the region enclosed by S. 



(C2) Positivity: 

Quasi-local energy is positive and is monotone in a suitable sense but vanishes in flat 
space-times. 



(C3) Limiting behaviors: 

(a) Standard sphere limit: Quasi- local energy recovers the standard value of mass en- 
closed in S for S ~ S 2 in a spherically symmetric space-time. In particular, for a sphere 
centered at the origin of Schwarzschild space-time, it coincides with the Schwarzschild mass 
parameter. 

(b) Marginally trapped surface limit: Quasi-local mass agrees with the irreducible mass 

/ Arca(5) 
V Wit ■ 
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(c) Small sphere limits: 

(i) in non-vacuum: on a small sphere of radius r centered at any arbitrary point p in 
the space-time M, the quasi-local energy-momentum recovers the energy-momentum of the 
matter observed by an equivalence class of instantaneous observers (the meaning of which is 
made precise in Sec. IIII CP characterized by a unit time-like and future-directed eo € T p M, 
namely, ^ L r 3 T(»,eo), where T is the energy-stress tensor. 

(ii) in vacuum: the quasi-local energy-momentum yields the analogue of the gravi- 
tational energy-momentum observed by the class of observers given in (i) in terms of the 
Bel- Robinson tensor [3] <B, namely, eo, eo, eo). 

(d) Large sphere limits: 

(i) at spatial infinity: quasi-local energy-momentum approaches the Arnowitt-Deser- 
Misner (ADM) energy-momentum in an asymptotically flat space-like hypersurface. 

(ii) at null infinity: quasi-local energy-momentum reproduces the standard Bondi mass 
Ebs and news flux J^-Ebs- 

In fact, in the large sphere limits, quasi-local quantities are no longer truly quasi- local as S 
contains an infinite measure. A perhaps more proper term here would be quasi- global. 

A similar scrutiny of the quasi-local angular momentum is technically more involved 
partially because the very definition of the quasi-local angular momentum in several contexts 
is yet to be unanimously agreed upon. Tentative investigations have been carried out in the 
past decades with few definitive outcomes (see 32) for an overview). The present work is 
thus focused on studying the quasi-local energy, only. 

Like any other construction of quasi-local energy, Liu-Yau's is subject to reasonable reality 
checks in order to be a physically sound candidate. Among (CI) -[C 3) listed above, (CI), 
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(C2), (C3)-(a), (b) and (d)-(i) have been discussed [2l|, [24j, 
other limiting behaviors of Liu-Yau quasi-local energy-the main body of the present work-is 
presented in the following sections. 

This article is structured as follows. In Sec. [Til a specific model of the quasi-local energy 

n 

proposed by Liu and Yau [21j is reviewed. The small- and large-sphere limits of the Liu- 
Yau quasi-local energy are closely examined, in Sec. [TTTJ and Sec. HVt respectively. Also 
considered, along the same lines, is the possibility of generalizing the notion of quasi-local 
energy for non-isolated gravitational sources. A summary is given in Sec. |Vj 
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II. DEFINITION OF THE LIU-YAU QUASI-LOCAL ENERGY 



The Liu-Yau quasi-local energy originated as a continuation of Yau's mathematical work 

36], sj]. The definition of the Liu-Yau quasi- local energy 



on the positivity of black hole mass 
is reviewed here for completeness. 

A. Physical part 

Consider a closed orientable space-like 2-surface S embedded in M. Vp G S, 3 null frame 
{Xtf^ G CT P M with its dual {0*}^ G CT p *M adapted to 5, such that CT P M = <CT P S © 
CTpS'- 1 , where CTpS 1 = span c {X 3 , Y 4 } and CTpS 1 - = span c {Xi, X 2 }. (Shorthand notation 
AE = E £g>]R A, where A is a field and E is a bundle, is used throughout. For example, 
for A = C, CE is the complexified bundle of E. When A = R, however, it is obvious 
that M,E = E.) The mean curvature vector of S at p in M is then H = — 2pX\ — 2pX2, 
where p = 9 l (D Xi X4) and p = 9 2 (D X4 X 3 ). It is worth noticing that the norm of H is, 
however, independent of the choice of moving frames. Indeed, = y/8pp, where pp > 
for space-like H. The physical part of the quasi-local energy is then defined as 



E phys (S) = — {\\H\\n 

87T J 



8tt 
s 



where Q is the volume form of S. 



B. Reference part 

Suppose that S, equipped with a Riemannian metric g s and Levi-Civita connection D s , 
has positive sectional curvature. Then by Weyl's embedding theorem |34 ] , there exists a 
unique isometric embedding i\ : (S,g s ,D s ) > (M°,g°,D°), up to isometries of R 3 , such 
that the second fundamental form 11° is solely determined by g and is positive definite on 
S. The composition of i\ with a successive embedding ii : R 3 Rf gives the reference 
embedding i = l 2 o i x : S ^ R 3 of S into Minkowski space-time M°. 

The same construction as in Sec. Ill Al gives the mean curvature vector, H°, of S at p 
in R 3 whose norm is \\H°\\ = ^8p°p°, where p°p° > 0. Hence the reference part of the 
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quasi-local energy is naturally defined as 

E rei (S) = — (\\H°\\Q. 
8n J 



A caveat is emphasized in [24] and [22[] to avoid any misleading interpretations. Unless S 
lies in a space-like hypersurface £ C M, it would be highly unnatural to require that S be 
iso metrically embedded into a space- like hypersurface in M°. That is, the absence of this 
additional hypothesis may result in positive quasi-local energy even in M°. 

One of the merits of the embedding scheme described above is that L\ is unique up to 
isometries of R 3 and that H° is therefore well-defined. However, it is in general non-trivial 
to obtain a complete solution to the full set of integrability conditions for the sequence 
of embeddings of S in M°. An alternative approach is to consider the co-dimension 2 
embedding l° : S M° at a possible expense of uniqueness (up to isometries of Rf) 
unless extra restrictions are imposed. The existence of t° is, nonetheless, guaranteed for any 
conformally flat S, as shown in 5|. In Sec. IIIIDt embeddings of this kind are realized as the 
null-cone reference. In Sec. IIVCI the asymptotic version of such embeddings is studied in 
the large-sphere limit at null infinity. 



C. The definition of the Liu-Yau quasi-local energy 
Definition 2.1 

The Liu-Yau quasi-local energy associated with the 2-surface S is 

E(S) =E rci (S) -E phys (S). 



III. THE SMALL-SPHERE LIMIT 



This section is devoted to gauging Liu-Yau quasi-local energy E(S) against criterion C3- 
(c) when S is a small sphere, as defined below. It is shown that E(S) satisfies C3-(c)-(i) for 
non-vacuum, but deviates, in vacuum, from the expected value in C3-(c)-(ii) by an extra 
term of which the physical nature is yet to be explored. 
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A. Construction of the small sphere 



The small sphere around an arbitrary point in the space-time is a space-like level set of the 
null cone emanating from that point. Vp G M, 3 a normal neighborhood U of p in M which 
uniquely determines a star-shaped neighborhood U of in T P M, such that the exponential 
map exp p is a diffeomorphism of U onto U whose inverse is denoted by exp" 1 . Without loss 
of generality, it is assumed that U is small enough such that the null cut locus C^(p) (£ U 
and hence that C^(p) (jL U. For a given orthonormal basis {ej}^ -1 for T P M (n — 4 when M 
represents a space-time), with (e^, ej) = S^ej, — 0, . . . , n — 1) and its dual basis {z i y~^~ l 
for T*M, the normal (Cartesian) co-ordinate system of the connection D on U is defined 
by £ = (x°, . . . , x n_1 ) G #(£/, R n ): x l = z l o exp -1 , i = 0, . . . , n — 1. Correspondingly, such 
a normal co-ordinate system (x°, . . . ,x n ~ 1 ) determined by {ej}]=o -1 assigns to each point 
q G U co-ordinates with respect to basis {e^}^) -1 of the pull-back exp" 1 ^) G U C T p M 

n-1 

via exp^ 1 (g) = x l (q)ei, Vg G £/. Thus, the normal (Cartesian) co-ordinate system 

i=0 

r -i i=n—l 

(x°, . . . ,x n_1 ) on [/ induces a Cartesian co-ordinate basis -I ^ > G 3^{U). 

Define the Lorentz radius function 5 G R) on M at p as 5(g) = |exp~ 1 (g)|, Vg G £/ 

f n— 1 

and consider the geodesic ball <5 _1 (c) = < q G £/ : Yl e i( xl (l)) = c 2 | C (7 in normal 

i=0 J 

co-ordinates for sufficiently small |c| > and a hypersurface (x°) _1 (t) for a given t G R. 
Then S(c,t) = <5 _1 (c) n (x°)~ 1 (t) = r _1 (Vc 2 - e t 2 ) is a closed sub-manifold of M, where 
r G [0, +oo)) by r = ^c 2 - e (x ) 2 is the radial co-ordinate of the spherical normal 
co-ordinates. In particular, for c = 0, 5 _1 (0) = J + (p,U) and thus Vg G U, 3 null Xi = 
exp~ 1 (g) G £/, such that the ray p : [0, a ) — > £/ by p(x°) = x°Xi, rr° G [0, ao), where 
[0, ao) is the maximal domain of p, uniquely determines a radial null geodesic from p to q, 
7 Xl : [0, a ) — >■ U by 7 Xl = exp p o p. The superscript X 1 stresses that 7 is the local flow of 
Xi, which is understood hereafter and the superscript will most often be suppressed when 
confusion is unlikely. It is noted in passing that as c = 0, x° = r, hence 7 is well affmely 
parameterized by r G [0, ao). The local null-cone A(p) = 5 _1 (0) — p is then foliated by 

= {7l(o,oo) : J™7( r ) = expp 1 ^).? e Hp)}- 
Definition 3.1 

For every given t, S(c,t) = S(0,r) = S 2 , as a regularly embedded space-like sub-manifold 
of M, is the desired small sphere. For economy of notation, S(0, r) will almost always be 
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abbreviated as S(r). 



B. Moving frames on the small sphere 

The regular embedding of S(r) into M suggests an envisaged choice of adapted or- 
thonormal null frames on U H A(p) as in Sec. [TTJ For every given r = x° G (0, ao), 
Vg G S{r) ^ M, 3 null frame {X<(g)}JzJ G CW 7 ,, where = T q M n T q M^ , with its 
dual {0*(g)}j=? G CH 7 ; adapted to 5, such that CT 9 M = CT q S(r) © CT 9 S'(r) ± , where 
CT q S(r) = span c {X 3 (g), X 4 (g)} and CT g 5'(r)- L = span c {Xi(g), X 2 (q)}. The principal 
bundle of null frames F(U fl A(p)) thus possesses a reduced structure group G = C x = 
GL(1,C) C SL(2,C), where the right group action F(U H A(p)) x G — ► F(C/ n A(p)) 
is given by u . z = (X[, X' 4 ) = (X x , . . . , X 4 )v4(z), where u = (X 1 , . . . , X 4 ) G CW 7 ,, 
Vg G U n A(p), z & G, and = diag(|^| 2 , 4j, |, |) G GL(1, C) = G. Such construction 
of adapted null frames is known, in physics literature, as the Geroch-Held-Penrose (GHP) 



formalism 14j, 11]. 



A few geometrical properties follow almost transparently from the preceding construction. 
For pedagogical purposes, however, it is considered helpful to first recall some facts about 
degenerate sub- manifolds of semi-Riemannian manifolds [20], which will be preliminary for 
Sec. [IV] as well. 

Denoted by AK the degenerate bundle over a degenerate sub-manifold (H, g^) of a semi- 
Riemannian manifold (M,g) equipped with a Levi-Civita connection D. It is known that 
AK = ATH n ATH 1 - and AK 1 - = ATH + AT if- 1 . In particular, when (H,g H ) is a null 
sub- manifold of a Lorentzian manifold (M,g), AK is the unique null A-line sub-bundle of 
ATH. Furthermore, if if is a null hypersurface, then AK = ATH- 1 . 
Definition 3.2 

(H, ph) is irrotational if DU G EndA ( AT H), VU G T(AK), where the bundle homomorphism 
DU G Hom A (ATii, ATM| H ), for a given U G T(AK), is defined as DU(X) = D X U, 
VX G ATH. 

The following properties are immediate consequences of the above series of definitions. 
Proposition 3.3 

(1) If (H,g n ) is irrotational, then DU G End A (ATii") is self-adjoint, VT7 G F(AK), i.e., 
gn(DU(X),Y) = g n (X,DU(Y)), VX, Y G ATH. 
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(2) If (H, g H ) is integrable and AK = ATH 1 - (or, equivalently, AK 1 - = ATH), then (H, g u ) 

is irrotational. 

Proof: 

(1) A straightforward calculation. 

(2) follows from (1). ■ 
Corollary 3.4 

Every degenerate hypersurface (H, gu) is irrotational. ■ 
Definition 3.5 

U G T(AK) is pregeodesic if 3/ G $(H, R), such that -D^f/ = fU. (H,gn) is geodesic if 
every [/ G r(Aif) is pregeodesic. 
Corollary 3.6 

If (if, <7h) is irrotational, then it is geodesic. 

Proof: compatibility and torsion-free properties of D. ■ 

In the present context, A(p) is a degenerate (in fact, null) hypersurface of M. Thus, with 
A = C, the following lemma, which will be used frequently throughout later calculations, is 
readily seen to hold. 
Lemma 3.7 

(1) e = ^(Dx.Xi) - i9 3 (D Xl X 3 ))' S ) = and k ee 6 2 (D Xl X 3 ) = 0. 

(2) p ee -^(D^Xg) and p ee ^ 1 ( j D X3 X 4 ) are real. 

(3) t - a - = and tt - a - = 0, where r ee -^ 2 (D X2 X 3 ), tt ee -6 2 (D X2 X 4 ) and 
57 + /3EE -6 2 (Dx 3 X 2 ) 

Proof: 

(1) Since Xi = 7 for the radial null geodesic 7, -D^Ax = 0. Hence the orthonormality of the 
null frame implies, for the Levi-Civita connection D, that 9 2 (Dx 1 X 2 ) = 0, 9 3 (Dx 1 X 3 ) = 0, 
and thus the claim. 

(2) Since A(p) is a null hypersurface of M, it follows immediately from Proposition 3.3 that 
^( J DA 1 (A 3 ))=^ 4 ( J DA 1 (A 4 )). 

(3) Similar to (2). ■ 

To better demonstrate the geometry of the local null-cone and the small sphere, it is 
instructive to introduce angular co-ordinate functions (1?, ip) G $(Un(A(p) — {70, 7^}), (0, tt) x 
[0, 27r)), where 70 and 7^ G ^\ are two disconnected leaves (aka generators) of J^a- Together 
with x° and r, the local null-cone (A(p) — {70,7^}) comes equipped with spherical normal 
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co-ordinates c = (x°, r, </?). Gauss' lemma then guarantees that |gf?r, £ <5T(?7n 

(A(p) — {7o,7tt})) constitutes a set of mutually orthogonal spherical co-ordinate basis. In 
particular, for every fixed r G (0,a ), the small sphere 5(r) can be parameterized by y?) 
or, as S* 2 = CP 1 , by the stereographic co-ordinates (CO) with the understanding that the 
poles are not covered. 

In the spherical normal co-ordinates <r, the null orthonormal frame at every point has co- 

n 

ordinate representation: Xi(q) = Yl X i(Q)^s , i = 1, • • • , n, V<? G U H (A(p) - {70, 7^}))- 



0=1 
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One noticeable simplification in the co-ordinate representation of Xi(q) is furnished by the 
foliation J^.. Recall from Sec. IIII Al that V7 G J^a, 7 = X 1; and that the affine parameter of 
7 can be arranged to be one of the spherical normal co-ordinates, namely q r = r. Therefore, 
along the local flow, r y Xl , of X\ in U, it is clear by definition that 1 = Dx ± r = X±r. Then, 
X* = 5® in the co-ordinate representation of Xi, and thus X 1 = 1 ■ 



C. Remarks on the vertex of the null cone 



Notice, by definition of normal co-ordinates, that x°(p) = r(p) = 0, whereas the angular 
part of q becomes degenerate at p. Consequently, the moving frames cannot be extended 
even continuously to p. It is then necessary to examine the directional dependence of the 
limiting process as r — ► 0. 
Definition 3.8 

The blow-up [l(| A(p) of A(p) at p is defined as A(p) = A(p) U n S 2 with the contraction 
7r : S 2 — > p. Hence A(p) and 5~ 1 (0) are homotopy equivalent. Correspondingly, the blow-up 
U of U at is simply the diffeomorphic pull-back of A(p) by exp" 1 , i.e., U = (U — 0) U^S" 2 , 
where 7r : S 2 — ► is the corresponding contraction in T p M. 

Now, for a given orthonormal basis {ej}*=Q _1 of T P M, it is natural to construct the 
spherical radial basis vector e r ($, (p) = sin$cos<£>ei + simfeiny^ + coshes, where ($,(p) G 
(0, 7i") x [0, 2n) represents the corresponding angular co-ordinates of S 2 obtained via the dual 
basis {z 1 }]^" of T*M. Then the direction-dependent limit of the null basis G 
T(CK) at S 2 coincides with X^p, 0) = e + e r , X 2 (p, 1?, 0) = |(e — e r ), X 3 (p) = ^(ei + 
ie 2 ), and X 4 = X 3 , V(tf, G S 2 . 

Physically, the blow-up of p defines an equivalence class of spatially isotropic instanta- 
neous observers at p. It suffices to consider only linear isometries of T P M, which constitute 
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a subgroup, Jz? = 0(3, 1, R), of the full isometry group 1(T P M) = if x R n . 
Definition 3.9 U 

An instantaneous observer at p is an ordered pair (p, eo), where eo € T p M is unit time-like 
and future-directed, (p, eo) is said to be spatially isotropic if the stabilizer of eo in if is 
if eo = 0(3,R). 
Definition 3.10 

Au equivalence class, denoted as ^ p (eo), of spatially isotropic instantaneous observers at p 
is defined as the orbit of if eo in T p M. 

It will be seen in Sec. IIII El that it is <^ p (eo) who measures the quasi-local energy in the 
small sphere limit at p. 

While S 2 is a topological 2-sphere, whether or not it can be realized in the small-sphere 
limit as a metric 2-sphere depends upon whether p is a curvature singularity, i.e., whether 
p G M. Analysis in the cases where p is of "elementary singularity" 26| is miserably 
complicated as no explicit Riemannian metric on the limiting sphere S 2 is available. Efforts 
have been channeled towards appealing to series expansions for small asphericity in S 2 , 
but have not proven to be as promising as expected. On the other hand, when the curvature 
at p is finite, it is always possible to choose S 2 as a metric 2-sphere as the limiting small 
sphere. Therefore, only the latter situation is considered in the present work. 



D. Reference embedding 

The reference part of the quasi-local energy associated with the 2-surface S(r) is solely 
determined by the intrinsic properties of S(r) although it appears to be defined in an extrinsic 
manner as in Sec. Ill Bl The apparent inconsistency is reconciled by virtue of the integrability 
conditions for S(r), particularly, as a sub-manifold of M. It turns out that for the current 
problem at hand, the Gauss equation alone is sufficient. 

As described in Sec. IIIBj consider the embedding i° : (S(r),g s ,D s ) (M°,g°,D°) of 
S(r) of dimension m = 2 into Minkowski (M°,g°) with flat connection D°. Now, \/p G 
S(r), with the usual identification of AT p S(r) and AT p 5'(r)- L with sub-spaces of AT t o p M, 
AT L o p M = AT p S(r) © AT p 5 , (r) _L . Here, it is assumed that T p S'(r)- L = span R {T, A^}, with 
orthonormal frame {T, N} and associated co-frame {8 T , 6 N }, where T is chosen to be time- 
like and N space-like. As a side remark, the standard space-time orthonormal frame field 
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is used in TM (A = R), whereas the null frame field is used in CTM (A = C). The same 
remark applies dually for the co-frames. The transformation between the frame {T, N} and 
the null frame {X 1: X 2 } as in Sec. IIIIBI is, by convention, fixed by (T,N) = (X 1 ,X 2 )il, 



where il 



1 1 

2 2 



1 -1 



G SO(l, 1, C). Together with any given orthonormal frame {-Ey}'! = i G 

AT p S(r) and its co-frame {o;- 7 }^"™ G AT*S(r), established are an adapted orthonormal 
frame {EjYjZi e AT L o p M and its co-frame {uj j } j pl G AT* 0p M satisfying 



E j 



and 



Alternatively, if null frames are used, 



Ea 



and 



1 


: 3 


i 

= 1 


N 


■ j 


= 2 


Ej 


■ j 


= 3,4 


9 T 


■ J 


= 1 


e N 


■ 3 


= 2 




■ 3 


= 3,4 


X j 


■ 3 


= 1,2 


X, 


■ 3 


= 3,4 


e j : 


3 


= 1,2 


& : 


3 


= 3,4 



Connections D s and D° are related by the shape tensor s -1 . By definition, \/x p ,y p G 
T p S(r), D°x p y = D h Xp y + s ± (x p ,y p ), where s ± (x p ,y p ) = II T (x p ,y p )T + II N (x p ,y p )N, with 
II T (x p ,y p ) = 9 T (D° p y) and II N (x p ,y p ) = 9 N (D° p y), where y is the extension of y p in the 
neighborhood of p. 

The fully contracted Gauss equation for the embedding (S(r),g s ,D s ) ^ (M°,g°,D°) 
can be written as 

4 4 

o = s s{r) + 6 2 [^^// iV (^, J E; J )/^(^)-(// Ar ) 2 



j=3 1=3 
3=3 1=3 



T\2 
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where e x = g{T, T), e 2 = g(N, N), and IF = £ LT (w', E t ), (r = N,T), or, for brevity, 

Z=3 



62 



:j/ iVx2 _ jjN . jjiV 



+ ei 



(JJ T ) 2 - II T ■ II T 



(1) 



Further simplification of the integrability conditions comes from the use of the so-called 
null-cone reference (aka light-cone reference). The same procedure as in Sec. IHI A[ applies 
for the construction of A°(p) in the reference space-time M° only now, it is a bona fide 
null cone as M° is Minkowski. One of the immediate benefits is that the leaves of ^ are 
shear-free, i.e., -a = 6 2 (D X ,X 3 ) = II N (X 3 , X 3 ) + II T (X 3 , X 3 ) = = -a° = 9 2 (D° Xi X 4 ) = 
II N (X 4 ,X 4 ) + II T (X 4 ,X 4 ). Then, Eq(CQ) becomes 



£2 



II N f - 1II N (X 3 , X 4 )U N (X 4 , X 



(II T ) 2 - 2II T (X 3 , X 4 )II T (X 4 , X 3 



■Ss(r) = 0. 

(2) 



Lemma 3.7-(2), when applied to A°(p), yields p° = p° and /i° = \i° . Hence, II N (X 3 , X 4 



-iff + 71° 



II N (X 4} X 3 ). Similarly, II T (X 3} X 4 



II (X 4 , X 3 ). In fact, finer results of p° and fi° can be obtained by basis transformation 
it p° = -6 2 (D X4 X 3 ) = -II T (X 4 , X 3 ) + II N (X 4 , X 3 ) = |( - II T + II N ) and, by the same 



token, fi° = \(II T +II N ). Hence, ||if || = y/8p^ = ^{II N ) 2 - (II T ) 2 . On the other hand, 
2II N (X 3 ,X 4 )II N (X 4 ,X 3 ) = \{II N f and 2II T {X 3 , X 4 )II T {X 4 , X 3 ) = \{II T f. Therefore, 
Eq(j2]) can be written as 



{ii N f - (ipy - 2s s{r) 



T\2 



which is the ultimate integrability condition needed here. It then follows that \\H°\ 



y/2Ss( r ), intrinsically determined by the (positive) sectional curvature of S(r) as desired. 



E. Limiting process 

As declared in Sec. IHI C\ it is always assumed that curvature at p is finite. Thus, the 
blow-up of p, S 2 , is a metric 2-sphere. Since the connection at p is flat, the leading term in 
every quantity is always its value in Minkowski space-time. The limiting process is then a 
straightforward but tedious radial expansion along the leaves of J^a of the quasi-local energy 
associated with the small sphere S(r) having a standard S* 2 as the limiting sphere at p. Such 
expansions have been carried out, up to various accuracy order, within the Newman-Penrose 
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(NP) formalism 25j and are well documented in the existing literature. In what follows, 
expansions of related quantities are quoted without proof as details can be found in, for 
example, Refs. 18], <J, fl. 1 

To differentiate among the behaviors of the quasi-local energy in non-vacuum and 
vacuum cases, the small-sphere limit is taken separately. 



1. Non- vacuum 

It turns out that the expansions of the relevant variables are needed only along the leaves 
of #a in the small-sphere limit. In the spherical normal co-ordinates <j on the local null 
cone of p, the following most pertinent quantities are expanded, in the Newman-Penrose 
formalism, in power series of the radial co-ordinate r. 



Ss(r) 
P 
/' 



2r- 2 + S^ + 0(r 



-r + -r> 00 



1 



S(r) 



-r 1 + 



tt = Vt Q r 2 



1 

-5 

3 



1 - -r 2 ^ + O(r 3 ) 



where f2 is the volume form of a metric 2-sphere Sq. Consequently, the quasi-local energy 
in non-vacuum can be written as 



E = ^ /(v 7 ^)- V8PP)V 

n S(r) 



r 
~2 



oo + So0?o + 3o A - 23o^io - 23 ^?) + 



So 

(Vo 2 + ^ + 2A° + ^S 1 -^ )+^ 



fio 

47T 



+ 0(r 4 ). 



Most of the terms in the integral vanish for one of two reasons shown below. The validity of 
the following lemma is well known, although it is sketched here with a slightly more formal 
proof. 



It is worth pointing out that Eq.(B18) in Q is a misprint, although the following Eqs.(B19) and (B20), 
are nevertheless, correct. 
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Lemma 3.11 

(a) J (3 /)fio = 0, V/ G tf*(-l) and J (9 #H = 0, Vg E ^ fc (+l), where ^ fc (s) is the sheaf 
So So 

of germs of spin weight s C-valued functions. 

(b) J Re^H = 0. 

So 

Proof: 



(1) An observation made in 12] indicates that, on So = CP 1 , the elliptic operator 9 is merely 



d in disguise. Hence, S : ^ fc (— 1) — > e tf k (0) is essentially the same as d : c^ 0,1 — > S 1,1 , 



where $ p,q is the sheaf of germs of C-valued forms of type (p, q) 33] . A similar argument 
holds for the conjugate operators. Then it is more transparent that the integral vanishes 
essentially by virtue of Stokes' theorem. 

(2) By definition, Re(V> 2 ) = -^(Wx^Xj, where W is the Weyl tensor. Then, using the 
basis transformation il at p, i.e., on S 2 , giWx^Xi, X 2 ) = ^^((eo + e r ) A (eo — e r ), (eo + 
e r ) A(e -e r )) = W(e Ae r ,e Ae r ) = a 2 W(e A e u e A e x ) + b 2 W(e A e 2 , e A e 2 ) + c 2 W(e A 
e 3 , e Ae 3 ) + 2a6# / (e Aei, e A e 2 ) + 2acy^(e Aei, e Ae 3 ) + 2bcW(e A e 2 , e Ae 3 ). Here, with 
a slight abuse of notation, W is also represented by W G SymBil(A 2 T p M x A 2 T p M, R), such 



that Kic(W) = 13]. In terms of the spherical harmonics YJ m ('i9, <£>), a = y^f C^i 1 — ^i 1 )) 
6 = {^.(yi+Yf 1 ), and c = y/^Y?, where G S 2 and {ejpj G T p M are as given in 

Sec. IIII CI Now it follows from the orthogonality of spherical harmonics and the trace-free 
character of the Weyl tensor that the integral vanishes identically. ■ 
Now, it is straightforward to establish the following: 

C fl , n _.n 2 ,„ 1 Vtn 



E{S) = r -\ [^o + 2A° + ±4,1 

So 

r 3 rl 1 
= y [-Ric(e ,e ) + -S 



J 47T 

+ 0(r 4 



° + 0(r 4 ) 



Air 

= _r 3 T(e ,e )+O(r 4 ), 

where T is the energy-stress tensor of matter. As anticipated, in non-vacuum, the leading 
contribution of E(S) in the small sphere limit at p comes from the energy of matter, ob- 
served by ^p(eo) (c.f. Definition 3.10). This is physically reasonable because any form of 
gravitational energy-momentum that is quadratic in curvature enters at higher orders in r. 
2. Vacuum 

In the vacuum case, higher order expansions are inevitably necessary. It is, nevertheless, 
a straightforward calculation to obtain the following expansions: 
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s. 



s{r) = 2r- 2 + Sfl ) +rSfl )+ r^l ) +0^) 
P = -r- 1 + ^r 3 ^° + 0(r 4 ) 



+ \r{^l + + \r 2 {^ 2 + ^2) + 



-0 



■0 



r U^Wo-77T g o(^i + # 1 )- 



,360 

J_~ / ,n— 

40 



40 



So(*i+4^)-^4M)+0(r 4 ), 



ft = tt r 2 



where 



S. 



S 



s 



(0) 
S(r) 

(1) 
S(r) 

(2) 
S(r) 



-5(80$! + Mi) 

-^(So^ + So^ 1 ) 
o 



In light of Lemma 3.11, the quasi-local energy in vacuum becomes 
E(S) 



1 

8^ 



S(r) 



S{r) 



72 
90 
90 



So 



So 
3 r 



<B(e , e , e , e ) - /(Re^ ) 2 ^ + 0(r 6 ) 



So 



05 (e , e , e , e ) - ^ (l£ 2 (ei, ei) + £ 2 (e u e 2 ) + (£ 2 (ei, e 3 ) + <£ 2 (e 2 , e 3 ) 



-€(e 2 ,e 2 )£(e3,e 3 )j+0(r 6 ), 

where 03 is the Bel-Robinson tensor [3j and (£(ej,ej) = — g(W eoej eo, ei), i,j = 1,2,3, is the 
symmetric (aka electric) part of the Weyl tensor. In the leading order, 0(r 5 ), the first term 
is the gravitational energy measured by p (eo), whereas the second term comes from the 
power series expansion in r of ^2Ss( T ) in the integrand of E(S). The physical interpretation 
of this additional term remains unclear. 
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IV. THE LARGE-SPHERE LIMIT AT NULL INFINITY 



For the sake of simplicity and also of highlighting the physics, the discussion is restricted 
to perfectly isolated sources in an empty (i.e., Ricci-flat) space-time. An attempt at gener- 
alization to non-isolated gravitational sources is briefly mentioned in Sec. IIVFI 

A. Construction of the large sphere at null infinity 

Several definitions and elementary properties of null infinity are collected here for the 
purpose of unifying terminology and notation. 

Definition 4.1 

A C r (r > 0) asymptote of a differentiable manifold M is an ordered triple (M,i,fo), 
where M is a manifold with boundary dM, / : M r — .\/ is embedding by inclusion, and 
fa : dM x [0, +oo) — ► M, such that f n (x, 0) = x, Vx G dM, is a collar [13]. 

The existence of collars in the differentiable category is easily shown [l7]. That boundaries 
of C° manifolds have collars is, however, far from obvious but is proved to be true pj. A 
collar fa can essentially be characterized, with recourse to the partition of unity, if necessary, 
by Q G $(U, [0, +oo)), having as its regular value so that fi -1 (0) = dM, where f/cMis 
a neighborhood of dM. Hence, a collar is denoted hereafter simply by Q. 

A richer structure of the asymptote becomes available when a differentiable manifold M 
possesses a semi-Riemannian metric g and a Levi-Civita connection D. In particular, in the 
Lorentzian category, it is a space-time, denoted by an ordered triple (M,g,D). The collar 
Q of the asymptote of a space-time then becomes crucially intertwined with the geometric 
structure. 
Definition 4.2 

A C r (r > 0) asymptote of a space-time (M, g, D) is a C r asymptote of M with a Lorentzian 
metric g and a Levi-Civita connection D associated to M. 

There exist in the literature various ways of defining the asymptotic structure of null 
infinity, most of which are essentially equivalent. It is attempted here yet another formulation 
that might be more appropriate in logic. Elaborate discussions of the topological properties 
of simple space-times can be found in, for example, , whereas more geometric properties 
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are recorded in, for example, 20] . For brevity, only future null infinity is considered for the 
past null infinity can be treated dually. 
Definition 4.3 

The null infinity (J?) of the space-time (M,g,D) is a (not necessarily connected) null sub- 
manifold of DM that is orientable and time-orientable. 

In fact, the time orientation of M, and thus of M, induces a compatible time orientation 
of J so that J = f+KSf-, where J !± = e ^n/ ± (M, M) = ^-I T (M, M) are, respectively, 
future {J^ + ) and past (^~) null infinity. Note that each of the J'^ = J" — J^ T is relatively 
clopen in J* and thus is a connected component of J? . 

In particular, J? has a collar, which, with a slight abuse of notation, is also denoted by 
Q, as null infinity is the only piece of DM that is of interest here. Moreover, the conformal 
properties of null structure in the Lorentzian category can be exploited to provide remarkable 
convenience in the analysis of the asymptotic structure of a space-time. 
Definition 4.4 

A space-time (M, g, D) is asymptotically empty and flat at null infinity if 

(1) It is conformally diffeomorphic to its asymptote at null infinity with i*g = Q 2 g and 
i*dQ ytz 0, where d is the exterior derivative in M; 

(2) 3 J^o = {q G J? : q is strongly causal} C J 1 ', which is thus open in M and is relatively 
open in J*\ 

(3) f2~ 2 Ric admits a C T extension to fi _1 (0); 

(4) JPq Pd S 2 x /, where / C R is a connected component of R. 

(M, g, D) is said to be asymptotically Minkowskian at null infinity when / = R. 

Definition 4.4 implies that J'q can be causally separated by a space-like surface S ~ S 2 C 
J^o- Hence there exists a tubular neighborhood jVi of S in J^o, i.e., there exists a line bundle 
(E, S, 7r s ), where E = TS"- 1 , and a C r diffeomorphism tj) : E — > jVi-, with ip(®x) = x, Wx G S, 
such that the diagram 

E 

commutes where ret : JVi — ► S is a retraction. Then, JV^ is foliated by curves 
that intersect S transversely. Furthermore, since J^o is a null hypersurface in M, it is 
irrotational and thus geodesic by Corollary 3.4 and 3.6. Hence, V future-directed null curve 
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A : $2 — > <^o, which can be parameterized to be a null geodesic, 3a G £2, such that A rh g S, 
where q = A(a). On the other hand, as A(a) is null, A(a) ^ T q S, which is space-like. Thus, 
either q £ S or q G S and A* a (T a <f 2 ) + T q S = T q J? . For notational purposes, denote X 2 = A, 
then X 2 G r(CL 2 ), where L 2 = TJ^ H T^- 1 is the only null line bundle over J^ - 111 fact, 
with a little hindsight, the normal bundle E over S might as well be chosen to be L 2 . 

Now, jV 2 can be co-ordinatized by g 2 = (u,$,ip) G #(^2,^2 x (0, 7r) x [0, 27r)), with the 
usual understanding that the poles of 5* 2 are not covered. Clearly S 2 is a non-empty locally 
acausal compact connected topological 2-sub-manifold of J?q, hence is sometimes called a 
"cut" of JV The set of all such cuts is denoted by <£y . 

To define the large sphere, first recall that a subset F C M is a causal (resp. chronological) 
future set in a space-time M if J + (F,M) C F (resp. I + (F,M) C F). For a given compact 
subset K C M, J + (K,M) is a causal future set of M. Consider A(if) = j + (K,M). It 
is shown in [23] that A(if) 7^ and that A(K) is a compact achronal embedded C° sub- 
manifold of M with boundary dA(K) = A(K) H o^ + 7^ 0- With the different iable structure 
prescribed above, A(K) is a closed null sub-manifold in M of co-dimension 1 with boundary 
dA(K). Therefore, A(K) has a tubular neighborhood jV\ in M. Similar to the treatment 
of JV 2 , A(K) is geodesic and JV\ is foliated by null geodesies 7 : — > A(K) such that 
X\ = 7 G r(CLi), where Li = 7M(iQ n IM(if) 1 is the only null line bundle over A(K). 
Definition 4.5 

The large sphere near future null infinity is a compact 2-surface K C M, such that A is a 
neat sub-manifold [l^] and that A (if) H J^q G 

According to Definition 4.5, Vs G cM, T X A(K) <f_ T x J"q ', i.e., A(if) is nowhere tangent 
to or, notationally, A(if) rtl J^ + . On the other hand, J° + has a collar, f2 say, which 
restricts to a collar, ^X(k)i 011 ^(-^O i n A(K). Hence, 3 a null foliation of JV\ each 
of whose generators, 7, is affinely parameterized by ^l^(K) so ^ na ^ = 1> where 

Xi = 7, and that i^i^ = £x 1 l P — 0, where £x 1 is the Lie derivative along the flow of X\. 
It is then natural to adopt g = (u,Q,$,<p) G $(<^o~,<o2 x [0, +00) x (0, tt) x [0, 27r)) as a 
co-ordinate chart on the neighborhood of in the asymptote. By the definition of a neat 
sub-manifold, A(K ) is covered by the chart (^, U) of M such that AnU = Q~~ l {u = uq), where 
uo G S 2 . Consequently, A admits an adapted co-ordinate system Qi = (uq, ^X(k)j ^> e 
S(I, [0,+oo) x (0,tt) x [0,2tt)). 

It is equally convenient to carry out the analysis in the asymptote using the local chart 
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g, as done in, for example, [23|. However, for the purposes of studying "physical fields" [151 ]. 
it is more often useful to work in the original space-time (M,g,D). The latter approach is 
taken in what follows. 

Most asymptotic behaviors of the original space-time near null infinity come almost for 
free because (M, g, D) is conformally diffeomorphic to its asymptote. For example, ^¥\V\M in 
M is foliated instead by null geodesies (possibly after reparameterization of null pregeodesics) 
7 : (6, +00) — ► M such that 7 = 207 j^j]. Let r be the affine parameter of 7, then 
Dx x r = 1. The local co-ordinate chart in M on the neighborhood of null infinity is simply 
q = (u,r,"&,ip), which is known, in physics literature, as Bondi-type co-ordinates, although 
in the original Bondi co-ordinates [4(, r is chosen to be a luminosity distance parameter as 
opposed to an affine parameter as used here. 



B. Moving frames on the large sphere 

The seemingly pedantic construction in Sec. IIV Al exhibits its advantages now when it 
comes to setting up adapted moving frames on A; almost the same moving frames as used in 
Sec. IIIIB1 can be applied in parallel for the large sphere K. The only difference lies in the 
obvious fact that rather than the blow-up sphere S 2 , the limit sphere now is K = Ad^^ ~ 
S 2 . Correspondingly, all quantities in the Newman-Penrose formalism are expanded in 
powers of -. 



C. Reference embedding 

Given a generic space-time (M,g,D), it is generally unlikely that a large sphere K in 
the asymptotic null region of M could be embedded into a genuine null cone in Minkowski 
reference space-time that is emanated from one single point. Thus the embedding scheme in 
Sec. IIII Dl becomes inappropriate in the large-sphere limit. However, it is possible to isomet- 
rically embed A into the asymptotic null region of the Minkowski space-time (M°,g°,D°) 
such that the sectional curvature of K, when calculated via the Gauss equation, is preserved 
regardless of which ambient manifold into which K is embedded. Next, recall from Sec. IIV Al 



that the leaves of the foliation of JV\ form a null congurence of X%. It is assumed 35J that the 



shear of such null congurence is the same at Kq in both M and M°. The same construction 
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as in Sec. IIVAI produces a local co-ordinate chart g° = (u, r°, ip) of the Bondi-type in M° 
on the neighborhood of J*§ that differs from g only by a possible reparameterization of the 
null congruence, registered by r° in g°. 

Asymptotic expansions of the Newman-Penrose variables in both M and M° can be 
performed at one stroke. In M°, ip° = 0, % = 0, ... ,4 for M° is flat. Some of the most 
pertinent expansions are listed below: 

-o 



<ra- (a°a ) 2 -Ua°^ + c.c.) 



P = ~ B \' TU ' — + 0(r 

_0 

°0 U , r,/ -3 



1 ^ + a^ + Qla° 

P = -7, 2 + °( r 

Zr r 

n = iv 2 (i-^-)+o ( r- 



.0 



1 n 2 7f°° 

where quantities with superscripts or subscripts represent their corresponding asymptotic 
values at Kq, which are not to be confused with those with superscripts o in the reference 
space-time (M°,g°,D°). 

Expansions of the full set of Newman-Penrose variables in both M and M°, when applied 
to the Gauss equation, establish an equality of the sectional curvature of K, calculated 
via two different embeddings, from which the relation between r and r° can be read off: 
r° = r + {d 2 a°° + %a° Q - 3go° - S^ ) + O^- 1 ). Now the assumption a°\ Ko = a o0 \ K o leads 
to a much simpler relation r° = r + 0(r _1 ). 



D. Bondi-mass loss 

After the preparatory work from previous sections, the calculation of the quasi-local 
energy in the large sphere limit now becomes completely transparent: 

E(K) = l-jl^/Sp^-^^ 

K 

= E BS + 0(r- 1 ), 
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where E-^s = — J (^2 + <T°a°)Qo is the Bondi mass loss. 
E. Energy flux 

The energy flux through K is defined as the rate of change in the quasi-local energy E(K) 
in the time-like direction characterized by T. Here, T(q) is related to {X i {q)} 1 0i G CTqK- 1 -, 
\/q e K by the same basis transformation il as in Sec. IIII Dl and agrees with the generator of 
the time translation subgroup of the Bondi-Metzner-Sachs (BMS) group at [28(. Hence, 



K 
K 

= -L j [y/S^pfj, - p) + 
((X 2 + ^X 1 ) V ^ 



87T J <9"U 
K 

in which t : K > M is the embedding map that induces t*(.£xi# 3 A 6 4 ) = —2p9 3 A 6 4 and 
l*(£x 2 3 A # 4 ) = 2/i6> 3 A 6> 4 by straightforward calculations. 
Very similarly, 



K 



8ir J du 

K 

where, again, r° = r + 0(r _1 ) is used in the co-ordinate representation of X\ and X 2 . 

Recall that the embedding scheme is tacitly chosen so that c 00 )^ = <7°\k and that 
■03 £ Then, with the help of the Newman-Penrose equations tp® + ~ ^°i J 4 — 

and ^4 + a° = 0, the flux of quasi-local energy in the large-sphere limit at J? + is 

dE(K) dE Tci (K) dE phys (K) 



dt dt dt 

1 

An 



Ko 

in which the leading term is precisely the flux of Bondi news. 
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F. Remarks on generalization to non-isolated systems 



Applying the notion of quasi-local quantities, in general, to non-isolated gravitational 
systems may incur failure to satisfy, for example, criterion C3-d. The obstruction to such 
generalizations largely lies in the difficulty of having a well-behaved or, well-described large- 
sphere limit. It suffices to analyze the problem at null infinity; the situation at the spatial 
infinity, if treated with care, is quite similar. Intuitively, it is conceivable that the non- 
isolated source has to eventually run off the boundariless manifold and wreck the topological 
structure of the asymptote described in Sec. IIV Al However, the following example illustrates 
a more subtle cause. 
Definition 4.6 

Given two quartic polynomials G(x) = 1 — x 2 — 2mAx 3 — e 2 A 2 x A and F(y) = —G(—y), 
where m > 0, e G R, and A > 0, assume that all of the roots of G(x) are distinct among 
which at least two are real (i.e., e / or mA < ^jp), denoted by x 2 < Xi, such that 
G(x) > 0,Vx G [x2,a;i]. The (charged) C-metrics, denoted by C\ = Ci(m, A,e), are a 3- 
parameter family of space-times of Petrov type-D, satisfying 

(1) topologically, C x = P\ x Qi, where P x = R x R+ and Q x 2* S 2 ; 

(2) metrically, (C\, gi, D\) is conformally diffeomorphic to its asymptote. In a local chart 
"^Ai = {t, U, x, z), set k to be a real constant, C\ ~ Pi x Q 1 , where (Pi, fi) is given by Pi « 
Rx[-x, +oo) and/i = -F{y)dt®dt+^0^, {Qxjii) is given by Q x rs (x 2 , Xi) x [0, 2ttk) and 
h x = ^f+G(x)dz®dz. Hence, g x = fi+h. The collar is defined by ft" 1 = A(x+y) G R+. 

The following lemma exposes one of the most peculiar features of the C-metrics, namely 
the so-called "conical singularity" or "nodal singularity" [ji3], [l| at one of the boundaries of 
the annulus Qi (in the proof, x\) that cannot be compactified in the C r category for r > 0. 
Lemma 4.7 

Q 1 2i S 2 only in the C° category but Qi = D 2 in the C r (r > 0) category, 
Proof: 

Consider another local chart vb x = (t, y, <p) which is C r -compatible with va 1: where 
r > 0, $0) = | -fe= and (p = k~ x z. Put p(tf) = y/G{x{$)) and $ = Then the 

graph of (i?, p(i?)) is (0, i? ) x (0,Pm] ; for some p M > 0. Clearly, p(0) = p($o) = 0. Now 
Qi~(0,tf )x[0,27r). 

(1) r = 0, simply extend continuously the domain of •& to compactify the annulus Q\. 
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(2) r > 0, the standard Bertrand-Puiseux test jscj, when applied on Q x , shows that [3] 
the C r -differentiable structure is preserved at the boundary of the annulus Qi if and only 



if k 



. However, 



loss of generality, set n 1 



dp 

dp 



unless m = or |e| = m > jj. Thus, without 





. Then, except for two 2-parameter families of electrovac 

i?=0 

solutions, m = or |e| = m > j-j, only one of the boundaries of the annulus can be 

compactified such that Qi = [x 2 , X\) x [0, 27r) = D 2 . ■ 

Lemma 4.7 clearly shows that a generic C- metric is not asymptotically empty and flat 

at null infinity in the sense of Definition 4.4. The remedy comes out of a key observation 

[l| that C\ is not maximal. Hence, a maximal extension of C\ leads to C = Ci U C2, 

where C 2 is an identical replicate of C\. In the asymptote of C 2 , C 2 ~ P 2 x Q 2 , where 

Lemma 4.7 applies except that in the local chart vb 2 of Q2 the opposite boundary of the 

annulus Q 2 is compactified, i.e., Q 2 — (x 2 ,Xi] x [0,27r) = D 2 in the C r category for r > 0. 

I I — 
Therefore, as shown in pj, in C, or its asymptote C, vb = {vb 1 ,v~b 2 } constitutes an atlas 

for Q = Qi U Q 2 = S 2 and that J" ~ S 2 x /, where / = R except for two null generators 

corresponding to precisely the boundaries of the annulus Q. But these two generators can 

be arranged such that they both are incomplete in the future or in the past. 

The above analysis demonstrates that for a non-isolated source modeled by the C-metrics, 

the topological and geometrical properties of J? as in Definition 4.4 can certainly be retained 

almost as for the isolated sources. The real difficulty, nonetheless, arises at a technical level 

when practical calculations are carried out. Evaluations of the Bondi mass or energy flux 

inevitably involve integrations on Q e which has to be covered by two charts. This task 

usually turns out to be analytically intractable [lo| . 



V. SUMMARY AND DISCUSSION 



The investigation of the limiting behavior of the Liu-Yau quasi-local energy is carried out. 
Such an analysis could be utilized to provide an appropriate certification for the Liu-Yau's 
proposal as a physically sound candidate for the quasi-local energy. Preliminary results that 
are considered new include: 

• In the small-sphere limit, the leading term in the quasi- local energy measured by the 
equivalence class of spatially-isotropic instantaneous observers ^ p (eo) at an arbitrary 
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point p in non- vacuum is found to be ^ L r 3 T(e , e ), where T is the energy-stress tensor 
of matter and r is the radius of the small sphere in the limit r —>■ 0. 

In vacuum, however, the gravitational quasi-local energy measured by & v gains 
an extra term in the leading order, in addition to the currently known limit 
^03 (e , e , e , e ), where 23 is the Bel- Robinson tensor. 

The occurrence of the extra term is considered as an example of how the quasi-local 
energy depends rather crucially upon the choice of the reference embedding. Since the 
co-dimension 2 embedding of the 2-surface S into the reference space-time M° is in 
general non-unique, it is plausible that an embedding scheme other than the null-cone 
reference may result in a different limiting behavior. Moreover, the currently known 
limit ^Q3(eo, eo, eo, eo) in vacuum is actually model dependent and usually variable 
to reference embedding (for example, Therefore, it is contemplated that the 

mismatch in the small-sphere limit in vacuum does not seem to serve as a strong piece 
of evidentiary support to rule out Liu-Yau's model. 

In the large-sphere limit at null infinity of an asymptotically empty and flat space- 
time, the Liu-Yau quasi-local energy is found to coincide, in radiating scenarios, with 
the Bondi mass loss and the news flux. A tentative generalization of the quasi-local 
energy to non-isolated sources encounters technical difficulties at null infinity in the 
example of the C-metric. 
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